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ABSTRACT 

A  method  of  numerical  solution  was  developed  for  obtaining  the 
stresses,  strains,  and  displacements  in  a  solid  propellant  grain  of 
finite  length  subject  to  loadings  of  pressure,  cure  and  thermal  shrink¬ 
age,  and  axial  acceleration.  Although  the  method  is  applicable  to  any 
solid  of  revolution,  the  development  was  restricted  to  a  propellant 
grain  in  the  form  of  hollow,  right,  circular  cylinder  bonded  at  the 
outer  periphery  to  a  rigid  motor  case.  The  ends  of  the  grain  were 
treated  for  both  bonded  and  free  conditions.  The  analysis  was 
limited  to  materials  which  are  homogeneous,  isotropic,  elastic,  and 
undergo  small  deformations.  Results  were  obtained  in  dimensionless 
form  which  can  be  expressed  in  terms  of  only  three  independent 
variables  -  Poisson' s  ratio  and  the  length-to-diameter  and  the  inner- 
to-outer  diameter  ratios  of  the  grain.  The  effect  of  these  variables 
on  the  stress,  strain,  and  displacement  distributions  in  the  pro¬ 
pellant  grain  was  examined. 

The  Raleigh-Ritz  method  was  employed  in  conjunction  with  the 
principle  of  minimum  potential  energy  to  obtain  expressions  for  the 
displacements  in  a  propellant  grain  in  the  form  of  a  hollow,  right, 
circular  cylinder  bonded  at  the  outer  periphery  to  a  rigid  motor  case 
and  subjected  to  cure  shrinkage.  The  analysis  was  limited  to 
materials  which  are  homogeneous,  isotropic,  elastic  and  undergo 
small  deformations.  Numerical  results  obtained  from  the  displace¬ 
ment  equations  compare  favorably  with  measured  values  and  those 
obtained  by  a  numerical  solution  using  relaxation  techniques. 
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DEFORMATIONS  AND  STRESSES 


IN  CASE-BONDED  SOLID  PROPELLANT  GRAINS  OF  FINITE  LENGTH 

BY  NUMERICAL  METHODS  ' 


Charles  H.  Parr 

Most  existing  stress-strain  analyses  of  solid  propellant  grains  have 
been  based  on  conditions  of  plane  strain,  which  is  equivalent  to  assuming 
the  grains  to  be  of  infinite  length.  These  solutions  are  applicable  to 
center  sections  of  very  long  cylinders  w/here  the  warping  of  cross 
sections  in  the  axial  direction  is  small.  These  plane  strain  solutions  for 
pressure,  thermal,  and  cure  shrinkage  loading  have  been  developed 
adequately  for  circular  cross-sections,  small  deformations,  linearly 
elastic  homogeneous  materials,  and  constant  material  properties  [  l]^. 
Significant  advances  have  been  made  in  plane  strain  analysis  to  include 
more  complex  cross  sectional  geometries  [2],  viscoelastic  properties 
[  3,  4] ,  and  large  deformations  [  4] . 

The  study  of  propellant  grains  of  finite  length  has  been  neglected 
until  recently.  Shortly  after  the  analysis  described  here  was  completed, 
two  papers  concerned  with  the  analysis  of  hollow,  right,  circular 
cylinders  using  numerical  methods  were  published  [  5,  6] .  Solutions  for 
stress  distributions  in  cylinders  of  finite  length  subject  to  various 
boundary  tractions  in  terms  of  continuous  functions  offer  formidable 
difficulties  even  in  classical  elasticity  theory  [  7,  8,  9j .  If  a  continuous 
solvition  can  be  obtained,  the  numerical  evaluation  of  the  solution  may  be 
quite  difficult,  particularly  in  certain  regions  of  the  cylinder. 

Following  several  abortive  attempts  at  this  Division  to  obtain  closed 
form  solutions  for  finite  length  grains  in  the  form  of  hollow,  right, 
circular  cylinders  with  the  required  boundary  conditions,  it  was  decided 
to  try  two  alternative  approaches:  (  1)  a  numerical  solution  of  the 
appropriate  differential  equations  and  (2)  approximate  solutions  based  on 
the  principles  of  minimum  energy.  This  section  of  this  report  is  based 
on  the  numerical  solution  method. 


Himbers  in  brackets  refer  to  references  at  the  end  of  the  report. 
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ANALYSIS 

Problem  Description 

The  most  general  end-effect  problem  in  the  mechanics  of  case- 
bonded  solid  propellant  grains  would  consider  a  system  of  three  concen¬ 
tric  cylinders  (motor  case,  liner,  and  propellant  grain) ,  the  inner 
cylinder  having  a  perforation  of  complicated  geometry  and  all  terminating 
in  a  complex  non- symmetric,  non- cylindrical  fashion.  The  end  condi¬ 
tions  would  vary  both  in  contour  and  in  the  conditions  of  restraint  on  the 
propellant  grain.  Loading  conditions  would  include  internal  pressuriza¬ 
tion,  cure  and  thermal  shrinkage,  three-dimensional  thermal  gradients, 
and  acceleration  forces  of  arbitrary  direction. 

In  addition,  the  propellant  under  some  circumstances  undergoes 
deformations  of  such  magnitude  that  the  theory  of  infinitesmial  deforma¬ 
tion  is  no  longer  valid  to  any  acceptable  degree  of  approximation;  the 
propellant  is  viscoelastic,  thus  making  the  problem  time  and  history 
dependent;  the  propellant  is  not  isotropic  in  its  response  to  stress;  and 
the  propellant  is  non-homogeneous. 

The  various  mathematical  theories  of  continuum  mechanics  are  at 

y 

such  a  stage  of  development  that,  with  the  exception  of  simultaneous 
large  deformations  and  viscoelastic  effects,  a  model  closely  describing 
the  actual  propellant  grain  can  be  formulated  in  differential  equation 
form.  However  the  complexity  of  the  complete  problem  i:f  such  that  ij:  is 
doubtful  that  existing  electronic  digital  computers  could  be  utilized  to 
obtain  numerical  results  in  reasonable  times.  More  serious,  however, 
is  the  lack  of  data  on  the  viscoelastic  response  of  propellants  under 
multiaxial  stress  conditions.  To  date,  there  is  not  a  single  propellant 
whose  viscoelastic  behavior  has  been  fully  characterized,  even  for  small 
deformations.  This  fact  alone  precludes  the  usefulness  of  the  more 
sophisticated  solution  techniques  available.  Lack  of  data  on  propellants 
dictates  a  choice  of  obtaining  a  complex  solution  in  terms  of  unmeasured 
physical  parameters  or  obtaining  a  simpler  solution  which  may  not  be 
compatible  with  reality. 

Accordingly,  to  reduce  this  problem  to  one  of  a  magnitude  solvable 
by  numerical  methods  on  existing  electronic  digital  computers  in  reason¬ 
able  times  and  wdiich  will  give  reasonable  qualitive  information  on  end 
effects  in  a  propellant  grain,  a  model  of  the  grain  in  the  form  of  a 
hollow,  right,  circular  cylinder  of  finite  length  with  the  outer  periphery 
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bonded  to  a  rigid  motor  case  and  with  the  ends  either  free  or  bonded  to 
the  case  was  considered.  The  cylinder  material  was  considered  to  be 
homogeneous,  isotropic,  and  elastic.  Small  deformation  theory  was 
used.  The  assumptions  of  case  rigidity,  i^iotropy,  and  small  deforma¬ 
tions  are  reasonable  if  (a)  the  interrelated  factors  of  interface  stress, 
case-propellant  modulus  ratio,  and  case  thickness  are  such  that  the 
deflections  of  the  motor  case  are  small  and  (b)  the  maximum  strains  in 
the  propellant  are  of  the  order  of  ten  percent  or  less.  These  conditions 
are  satisfied,  in  general,  by  motors  which  undergo  thermal  or  cure  ^ 
shrinkage  and  acceleration  loads  but,  in  general,  are  not  met  by  motors 
undergoing  internal  pressurization,  in  which  case  the  validity  of  these 
assumption  must  be  established  in  individual  problems.  This  can  be 
done  by  studying  the  plane  strain  solutions  for  both  elastic  and  rigid 
motor  cases.  The  assumption  of  homogeneity  is  adequate  if  stresses 
are  only  considered  on  the  macroscopic  level.  The  assumption  of 
elasticity  was  made  for  simplicity  although  the  extension  of  this  method 
to  include  viscoelastic  effects  is  under  study. 

Problem  Formulation 

The  model  used  in  this  analysis  is  shown  in  Fig,  1,  It  consists  of  a 
hollow,  right,  circular  cylinder  of  finite  length  which  is  bonded  at  the 
outer  periphery  to  a  rigid  solid.  The  cylinder  material  was  homogeneous, 
elastic,  and  isotropic.  In  the  present  study  two  types  of  end  conditions 
were  considered:  (a)  both  ends  free  and  (b)  one  end  bonded  to  a  rigid 
wall  and  the  other  end  free.  For  cylinders  with  free  ends  under  thermal 
or  cure  shrinkage  loading  or  pressure  loading  an  additional  symmetry  is 
introduced  so  that  the  numerical  solution  was  carried  out  only  for  one- 
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half  of  the  length. 

The  cylinder  was  considered  to  be  loaded  by.internal  pressure  on 
the  inner  radius  and  on  each  of  the  free  ends;  by  a  linear  expansion,  6, 
which  would  havg  a  negative  value  for  shrinkage;  and  axial  acceleration. 

Stress  functions  introduced  by  Southwell  [  10],  but  which  are 
described  more  clearly  by  Allen  [  ll],  were  used  for  the  numerical  solu¬ 
tion.  These  Southwell  stress  functions  are  particularly  well  adapted  to 
the  use  of  relaxation  methods  for  the  numerical  solution  of  problems 
involving  axial  symmetry.  Southwell  shows  that  if  two  functions  r,  z) 
and:  '*'(  r,  si)  satisfy  the  second  order  differential  equations 
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the  compatability  relations  and  equilibrium  equation  are  satisfied  for  an 
elastic  solid  undergoing  small  deformations  if  the  stresses  are  defined  by 
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where  r,  0,  z  denote  a  cylindrical  polar  coordinate  system,  a  denotes 
stress,  and  v  denotes  Poisson*  s  ratio.  Eqs.  1  through  6  were  general¬ 
ized  to  include  the  effects  of  body  forces  due  to  axial  acceleration  and 
propellant  expansion  due  to  temperature  change  and  cure  shrinkage. 

The  addition  of  the  expansion  term  has  been  discussed  previously  by  Parr 
[12].  Acceleration  terms  are  discussed  by  Southwell  [  13] .  If  the  pro¬ 
pellant  cylinder  is  under  an  axial  acceleration  A  and  undergoes  a  linear 
expansion  6,  the  stresses  and  strains  may  still  be  defined  in  terms  of 
stress  functions  which  satisfy  Eqs.  1  and  2.  For  convenience,  Eqs.  1 
through  6  and  the  allied  strain  and  displacement  equations  were  put  into 
dimensionless  form.  Defining 
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where  u  indicates  displacement,  €  indicates  strain,  b  is  the  outer  radius 
of  the  cylinder,  y  is  the  density,  and  E  is  the  elastic  modulus  of  the 
cylinder,  the  partial  differential  equations  for  the  stress  function  p,  tj) 
and  p,  nl  may  be  written 

-  1  +  1-1  r  0 
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The  displacements,  strains,  and  stresses  are  defined  by 
u  =  ii-1-y.l  t|j  +  [l-v](l>j  +pf6-  vpf)Z  , 
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The  function  f{  p)  appearing  in  Eq.  10  will  be  discussed  later.  The  dimen¬ 
sionless  coordinate  system  is  illustrated  in  Fig.  2.  The  origin  of  the  r\ 
coordinate  is  taken  at  the  motor  midsection  for  symmetrical  problems 
and  at  the  left  end  of  the  cylinder  for  unsymmetrical  problems. 

Boundary  Conditions 

In  addition  to  satisfying  Eqs.  7  and  8,  the  stress  functions  (}>  and 
have  certain  restrictions  placed  on  them  at  the  boundaries  which  are 
determined  from  the  boundary  conditions. 

Boundary  p  s  o* 

On  the  inner  boundary  of  the  cylinder,  r  =  a,  the  most  general  form 
of  boundary  conditions  for  the  problem  under  consideration  can  be 
expressed  as 

=  -  Pzl  z) 

and 

where  p2(  z)  is  the  normal  pressure  on  the  boundary  and  t^(  z)  is  the 
shearing  stress  on  the  boundary.  In  dimensionless  form,  the  inner 
boundary  is  ^  ;  a  s  ^  and  the  boundary  conditions  become 

(n) 

and 

Boundary  n  :  H 

The  free  end  of  the  cylinder  was  considered  to  be  loaded  by  normal 
pressure  p3(  r)  and  shear  stress  t3(  r) .  In  dimensionless  form  this  * 
boundary  becomes  tj  r  H  and  the  boundary  conditions  can  be  expressed  as 

%  =  =  -  Hslp) 

and 

Boundary  p  z  1 

On  the  boundary  r  =  b,  or  p  r  1  in  dimensionless  form,  the  bonding  o 
to  the  rigid  motor  case  requires  zero  radial  and  axial  displacements. 

Hence  in  dimensionless  form  these  boundary  conditions  become 
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1.  Symmetrical  Problems 


F or  symmetrical  problems  the  plane  «  0  was  taken  as  the  mid¬ 
section  of  the  grain.  From  considerations  of  symmetry  on  the  stresses 
and  displacements  it  can  be  deduced  that  the  symmetry  boundary  condi¬ 
tions  are 


2.  Unsymmetrical  Problems 

^  For  unsymmetrical  problems  the  plane  q  :  0  was  taken  at  the  left 
end  of  the  cylinder. 

(a)  When  this  end  is  bonded  to  a  rigid  case  the  radial,  and  axial 
displacements  are 

u  r  0 
P 

and 

( b)  When  this  end  is  free,  the  boundary  conditions  become 
analogous  to  Eqs.  23  and  24,  i.e., 

*  Hi  (  p) 

and 


»p,=  T.{p)  . 

If  the  axial  displacement  is  zero  along  the  coordinate  t|  x  0,  the 
function  f(  p)  in  Eq.  10  is  zero.  If  the  axial  displacement  is  not  zero 
along  this  coordinate,  the  function  f(p)  must  be  separately  evaluated. 
Since  such  cases  are  not  considered  in  the  results  covered  by  this  report, 
the  evaluation  of  f(  p)  will  not  be  discussed  further. 

For  the  loadings  considered  in  this  report,  the  normal  surface  loads 
n.  are  considered  constant  and  the  shear  loads  T.  are  all  taken  to  be  zero 
The  extension  of  the  finite  difference  forms  to  variable  loads  is  quite 
B  t  raightf 0  rw'a  rd,  however. 

The  reduction  of  Eqs.  7  through  32  to  finite  difference  form  and  the 
method  of  solution  is  outlined  in  the  Appendix. 
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RESULTS  OF  COMPUTATIONS 

An  inspection  of  the  field  equations  of  the  stress  function  (Eqs.  7 
and  8) ,  the  definition  of  dimensionless  displacements,  strains,  and 
stresses  (Eqs.  9  through  18),  and  the  boundary  conditions  (Eqs.  21 
through  32)  reveals  that  the  dimensionless  displacements,  strains,  and 
stresses  are  linear  functions  of  the  load  parameters  6,  Z,  and  T.  , 
Thus  the  parameters  whose  effects  need  investigation  under  each  type  of 
loading  are  Poisson' s  ratio,  v,  and  two  geometric  variables  which  were 
somewhat  arbitrarily  chosen  to  be  the  ratio  of  cylinder  length-to-diameter, 
\,  and  the  outer-to-inner  radius  ratio,  K. 

Some  90  problems  with  variations  of  these  three  parameters  for 
various  loads  have  been  worked  numerically.  Tabulation  of  the  variation 
of  the  four  stresses  and  two  displacements  in  critical  areas  as  a  function 
of  these  three  parameters  would  require  not  only  a  discourse  of  consider¬ 
able  proportions  but  a  massive  amount  of  data  reduction.  Accordingly, 
the  presentation  of  the  results  of  these  calculations  is  being  prepared  in 
two  parts.  This  report  will  primarily  be  concerned  with  overall  deforma¬ 
tion  effects  and  stresses  at  the  midsection  of  such  motors.  This  is  in 
line  with  the  interest  of  this  Division  in  profile  measurements  as  a  non¬ 
destructive  test  and  in  the  use  of  motors  with  propellant  grains  of  the  same 
geometry  employed  in  this  analysis  for  the  determination  of  failure 
criteria. 

The  second  part  of  the  presentation  of  results  consists  of  a  more  com¬ 
plete  stress  study  and  is  being  prepared  by  John  Wise  of  Thiokol  Chemical 
Corporation,  Redstone  Division,  in  the  form  of  design  charts  suitable  for 
use  by  rocket  motor  designers.  This  study  will  be  published  separately. 
Shrinkage  Loading 

Fig.  3  shows  the  theoretical  radial  displacement  of  the  inner  radius 
for  a  grain  with  free  ends  33  in.  in  length  having  an  outer  radius  of  3  in. 
and  an  inner  radius  of  1  in.  for  four  sets  of  Poisson' s  ratio  and  linear 
expansion  factor.  Superposed  are  measured  displacements  due  to  ther¬ 
mal  and  cure  shrinkage  from  a  6C2-33  static  test  motor  (as  1,  b  s  3, 

I  s  33) ,  made  with  nitrocellulose  plastisoi  propellant.  The 

measurements  were  made  of  the  profile  of  the  propellant  in  the  static  test 
motor  had  ceased  to  change  with  time,  indicating  that  cure  shrinkage  was 
complete  and  viscoelastic  effects  were  no  longer  affecting  the  profile. 
Neither  of  the  properties  v  or  6  is  known  accurately  for  this  propellant. 

it' 
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The  computed  displacements  are  directly  proportional  to  6,  and  therefore 
for  a  given  v  it  is  possible  to  select  a  value  of  6  which  fits  the  experi¬ 
mental  data  at  the  midsection  of  the  motor. 

Accordingly,  displacement  curves  for  v  a  0.49  and  v  a  0.50  were 
fitted  to  the  actual  displacements  at  the  midsection  of  the  motor,  giving 
values  of  6  of  -  0. 0050  and  -  0. 0043,  respectively.  ( For  this  motor, 
the  radial  displacement  at  the  center  for  v  a  0.5  is  98.  5%  of  the  plane 
strain  solution.)  To  demonstrate  the  effects  of  variations  in  Poisson' s 
ratio  and  the  linear  expansion  factor  6,  curves  are  also  shown  for 
V  s  0.  5,  6  ■  -  0. 0050  and  v  a  0. 49,  6  s  -  0. 0043. 

It  can  be  seen  that  any  value  of  Poisson' s  ratio  from  0.  50  to  0.49 
(and  perhaps  lower)  when  combined  with  a  suitable  value  of  6  will  give  a 
reasonable  fit  to  the  experimental  data.  This  suggests  a  method  of  deter¬ 
mining  Poisson' s  ratio  if  the  unconstrained  value  of  the  shrinkage  is 
known.  The  converse,  however,  would  not  apply  since  small  errors  in 
the  value  of  Poisson' s  ratio  would  lead  to  large  errors  in  the  shrinkage. 
This  demonstrates  the  dependence  of  the  apparent  shrinkage  on  the 
dilatancy  or  compressibility  of  the  material. 

Figs.  4  and  5  demonstrate  the  distortion  produced  in  the  longitudinal 
cross  section  of  a  long  and  a  short  grain  under  shrinkage  loading.  These 
grains  had  both  ends  free  so  that  the  deformations  were  symmetrical 
about  the  transverse  centroidal  plane.  Fig.  6  demonstrates  the  distortion 
in  a  propellant  grain  with  the  left  end  fixed  and  Fig,  7  shows  the  corres¬ 
ponding  dimensionless  contour  map  of  the  stresses.  It  can  be  seen  from 
Figs.  6  and  7  that  there  is  little  interaction  of  end  effects  so  that  the  left 
side  may  be  taken  as  being  typical  of  rigidly  bonded  ends  and  the  right 
side  as  being  typical  of  free  ends. 

Figs.  8  through  10  demonstrate,  respectively,  the  effect  on  the  pro¬ 
file  of  Poisson' s  ratio,  v,  length-to-diameter  ratio.  A,  and  radius  ratio, 
fC,  on  grains  with  free  ends  subjected  to  shrinkage  loading. 

The  examples  in  Figs.  4  through  10  show  results  of  the  type  that  can 
be  obtained  with  the  solution  method  outlined  here.  A  more  detailed 
study  was  made  to  determine  the  values  of  Aj  the  length-to-diameter 
ratio  which  was  needed  to  produce  a  condition  of  plane  strain  at  the  mid¬ 
section  of  the  propellant  grain  with  free  ends.  Figs.  11  through  15 
demonstrate  the  effect  of  the  parameters  A,  K,  and  v  on  the  deviation  of 
the  radial  displacement  of  the  inside  surface  at  the  middle  of  various 
cylinders  from  the  displacement  of  the  inside  surface  in  a  cylinder  of 


l! 
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infinite  length.  The  dashed  lines  indicate  areas  where  the  data  was 
insufficient  to  give  reliable  values.  The  ratio  of  the  radial  displacement 
at  the  inside  radius  to  the  inside  radius,  ~  ,  in  a  cylinder  of  infinite 
length  is  given  [15]  by 

\  [l  +  v][/<^-l]6  ^  (33) 

^  '  l+[l-2v]/<* 

From  Figs.  11  through  15,  it  can  be  seen  that  the  end  effects  are 
most  extensive  for  v  r  0.5  and  decrease  rapidly  with  decreasing  v.  End 
effects  become  less  extensive  with  decreasing  radius  ratio  K  (decreasing 
web  fraction) .  Since  the  end  effects  were  most  extensive  for  v  s  0.  5,  the 
largest  number  of  calculations  were  made  for  this  value  of  Poisson' s 

ratio. 

By  cross  plotting  Figs.  11  through  15,  the  ratio  of  radial  displace¬ 
ments  can  be  plotted  as  a  function  of  reduced  web  w’'''  s  (b  -  a)  /2b,  which 
is  sometimes  a  more  convenient  parameter  to  use.  A  typical  plot  of  this 
type  is  shown  in  Fig.  l6  for  v  =  0.5  for  various  values  of  The  curve  ^ 

\  =  0  in  Fig.  16  is  plotted  from  the  theoretical  plane  stress  solution.  Note 
that  it  does  not  reach  '100%  of  the  plane  strain  solution  at  w*  =  0.  The 
assumptions  of  both  plane  strain  and  zero  length  are  untenable  and  lead  ^ 
to  the  discontinuity  shown  here. 

Figs.  17  through  22  show  the  circumferential  and  axial  stresses  at 
the  midsection  of  the  grain  with  free  ends  as  a  function  of  v ,  A.,  and  K. 

The  dashed  lines  indicate  areas  where  the  data  was  insufficient  to  give 
reliable  values  of  the  stresses.  The  compressive  values  of  axial  stress 
for  \  <  1  in  Fig.  18,  while  surprising,  are  apparently  correct.  Several 
points  were  calculated  in  this  area  to  establish  this  behavior.  •  ^ 

One  proposed  method  of  testing  propellaifts  in  a  multiaxial  stress  •  ^ 

field  utilizes  a  rocket  motor  similar.in  all  respects  to  the  model  dis¬ 
cussed,  with  loading  introduced  by  lowering  of  the  temperature  of  the 
system.  The  differential  thermal  expansion  of  the  case  and  propellant 
results  in  the  type  of  shrinkage  considered  here.  It  is  bf  interest  to  note 
since  Poisson' s  ratio,  v,  is  a  function  of  temperature,  t|;^at  as  the 
temperature  changes,  not  only  will  the  displacerhent  in  a  motor  o|  given 
dimensions  change  but  the  shape  of  the  displacement  curve  w^l  also 
change,  indicating  that  the  stress  field  is  also  changing  in  form.  Ihus 
while  Figs.  18  and  19  for  v  :  0.  5  may  be  valid  Jo  r  a*  failure  at  room 
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temperature,  the  stresses  at  failure  at  a  temperature  of  -  40“F  ne«id  to 
be  determined  from  another  curve,  say  for  v  z  0.40. 

Acceleration  Loading 

To  date  an  extensive  study  of  acceleration  loads  has  not  been  made; 
however,  such  a  study  is  in  progress.  Typical  distortion  plots  are 
shown  in  Fig.  23  for  a  grain  with  both  ends  free  and  in  Fig.  24  for  a 
grain  with  one  end  free  and  one  end  fixed.  The  stress  contours  for  the 
grain  of  Fig.  24  are  s^iown  in  Fig.  25.  Again,  the  end  effects  are  well 
separated  so  that  the  right  end  is  typical  of  a  free  end  and  the  left  end  is 
typical  of  a  fixed  end. 

One  striking  property  of  the  cylinder  with  both  ends  free  is  that  the 
stresses  and  defections  are  antisymmetric  about  the  midsection.  An 
interesting  result  of  this  fact  is  that  the  stresses  and  deflections  at  the 
midsection  of  any  axially  accelerated  cylinder  of  finite  length  are  exactly 
the  same  as  the  stresses  and  displacements  everywhere  in  a  similar 
cylinder  of  infinite  length. 

Pressure  Loading 

A  few  cases  of  pressure  loading  have  been  run.  However,  the  use¬ 
fulness  of  this  solution  is  limited  to  relatively  thick  and  stiff  motor  cases 
since  the  model  assumes  a  rigid  case.  An  additional  limitation  is 
imposed  since  the  magnitude  of  strains  under  usual  operation  pressure 
will  sometimes  exceed  the  range  of  validity  of  the  analysis.  Details  of 
the  investigation  for  this  loading  will  be  presented  at  a  later  date. 

DISCUSSION 

Shrinkage  Loading 

It  was  determined  that  within  the  range  of  allowable  strains  in  which 
the  analysis  is  valid,  say  10%  or  less,  the  end  effects,  expressed  in 
dimensionless  form  are  functions  of  only  Poisson's  ratio,  v  ,  the  ratio  of 
outer-to-inner  radius,  K,  and  the  length-to-diameter  ratio,  All  three 
parameters  affect  the  nature  of  the  end  effects  significantly  in  their 
ranges  of  interest. 

Only  a  few  profiles  of  propellant  grains  are  presented  in  this  report; 
however,  comparison  of  the  many  profiles  which  have  been  calculated 
showed  that  the  apparent  flatness  of  the  deflection  curve  in  the  center 
section  of  the  cylinder  is  not  a  reliable  indication  of  the  degree  of 
approach  to  the  plane  strain  condition.  In  actual  measurements,  the 
apparent  flatness  is  further  clouded  by  normal  scatter  of  data. 
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Consideration  of  these  factors,  indicates  that  the  use  of  case-bonded 
propellant  cylinders  for  obtaining  values  of  the  thermal  coefficient  of 
expansion  and  of  stress  and  strain  at  rupture  is  not  at  all  straightforward 
and  may  lead  to  large  errors  if  a  plane  strain  analysis  is  used.  For 
instance,  suppose  a  motor  with  it  =  6  and  \  s  4  is  used  to  predict  the  ther¬ 
mal  coefficient  of  expansion.  For  v  =  0.  50  the  deflection  expected  can  be 
seen  from  Fig.  3  to  be  only  65%  of  the  plane  strain  deflection.  If  the  pre¬ 
diction  is  based  on  a  plane  strain  solution,  the  apparent  coefficient  of 
thermal  expansion  would  only  be  65%  of  the  true  coefficient  of  thermal 
expansion.  If,  however,  v  is  known  to  be  0.49,  the  apparent  coefficient 
of  thermal  expansion  would  be  84%  of  the  true  value.  However,  if  v  is 
actually  0. 49  for  the  propellant  but  it  is  assumed  to  be  0.  50,  the  apparent 
coefficient  of  thermal  expansion  using  Eq.  33  would  be 

(Plane  strain  displacement  for  v  « 

84%  X  I'Piane  strain  displacement  for  v  a  O.ouj 

or  48.  5%  of  the  true  value.  Errors  of  a  similar  type  could  occur  in  the 
prediction  of  stresses  and  axial  strains. 

Viscoelastic  Effects 

The  elastic  solutions  obtained  here  are  applicable  to  viscoelastic 
materials  at  long  times  if  the  equilibrium  moduli  are  used  or  for  very 
short  times  if  the  glassy  moduli  are  used.  Since  the  tensile  modulus 
occurs  only  in  the  dimensionless  stress,  pressure,  and  acceleration 
parameters,  more  general  conclusions  can  be  reached  in  certain  cases. 
For  instance  with  time  independent  Poisson' s  ratio,  the  viscoelastic 
strain  and  displacement  solutions  for  shrinkage  loading  are  the  same  as 
the  elastic  solutions.  Further,  although  a,  6,  and  E  are  functions  of 
time,  the  dimensionless  stresses  (ct/-6E)  are  constant,  allowing  the 
stresses  to  be  evaluated  if  the  time  functionality  of  6  and  E  are  known. 

Additional  viscoelastic  solutions  can  be  obtained  by  using  the  calcu¬ 
lated  numerical  results  for  the  elastic  solution  and  utilizing  numerical 
methods  or  other  approximate  methods  for  inversion  of  Laplace  trans¬ 
forms.  Such  methods  are  under  study  and  results  will  be  reported  at  a 

later  date. 

Curved  Ends 

A  computer  program  has  been  completed  to  investigate  grains  with 
curved  ends  which  are  bonded  to  the  motor  case.  It  is  expected  that  the 
effects  of  the  shape  of  the  curved  end  on  propellant  stresses,  strains, 
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and  displacements  will  be  studied  using  this  program.  This  will  be 

reported  at  a  later  date. 
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APPENPIX 


SOLUTION  BY  RELAXATION  METHODS 


Finite  Difference  Approximations 

For  the  numerical  solution  of  Eqs.  7  and  8,  subject  to  appropiriate 
boundary  conditions,  the  cylinder  cross  section  in  the  p,  i)  plane  was 
divided  into  a  square  mesh  as  illustrated  in  Fig.  2.  Each  intersection  of 
the  mesh  was  used  as  a  nodal  point  for  finite  difference  formulation. 

Using  three-point  central  difference  approximations  from  Taylor' s  series 
[  11],  the  following  approximations  result. 


and 

^  . 

8p*  Ap  * 


(A-1) 

(A-2) 


Similar  expressions  can  be  derived  involving  ij'  andq  ,  With  the  square 
mesh  used,  Ap  =  Aq=  h. 

Numeric  subscripts  corresponding  to  Fig.  26  were  used  to  denote 
node  location  with  respect  to  the  central  node.  It  should  be  understood 
that  this  notation  applies  to  the  neighborhood  of  the  node  of  interest  and 
is  independent  of  the  notation  that  locates  a  particular  node  in  the  body, 
Applying  relations  of  the  types  Eq.  A-1  and  A-2,  Eqs.  7  and  8  can 
be  expressed  in  finite  difference  form  as 


[(|>4  -  4*2]  +  (|>1  +  4>2  +  4*3  +  ‘t'i  “  'M'o  =  0  (A-3) 

and 

h 

77  [  “  4^2]  +  4^1  +  4^2  +  4^3  +  4^4  -  44^0  -  4>i  ”  4>3  +  24>o  :  0  .  { A-4) 

The  values  of  the  stress  functions  at  points  in  the  interior  of  the  body 
must  satisfy  the  field  equations  in  finite  difference  form,  Eq.  A-3  and  A-4. 

Finite  Difference  Form  of  the  Boundary  Conditions 

In  writing  the  central  finite  difference  approximation  for  a  normal 
gradient  at  a  boundary,  the  resulting  expression  involves  nodal  points 
lying  outside  the  boundary  which  are  called  fictitious  nodes.  By 
simultaneous  application  of  the  finite  difference  forms  for  the  governing 
field  equations  and  the  boundary  condition  equations  at  the  nodes  on  the 
boundary,  these  fictitious  nodes  can  be  eliminated.  This  is  the  standard 
procedure  for  hand  relaxation  as  used  by  Southwell  and  others  [  11,  13] . 


Boundary  p  s  Of 

The  boundary  conditions  to  be  considered  are  Eqs.  21  and  22  where 

jSi  is  a  constant  and  Ta  :  0  .  Thus 

a  5  -  n2 
P 

and 


Expressing  these  stresses  in  terms  of  the  finite  difference  forms  by  the 
use  of  EqS,  15  and  18  and  relations  of  type  expressed  by  Eq.  A-1,  the 
boundary  conditions  Eqs.  A-5  and  A-6  become 

<l>2  -  <|»4  +  >1^2  -  4*4  -  2Q'hn2  =  0 

and 


i|ij  -  (jjl  s  0  . 

However  since  the  relation 

84 


IS 


:  0 


resulting  from'  Eq.  18  and  £q.  A.6  is  valid  at  constant  p,  it  may  be  ■ 
integrated  directly  along  the  boundary  to  yield 
:  constant  . 

Since  there  is  one  degree  of  arbitrariness  in  the  stress  functions  (see 
Allen  [  1 1]  I  p.  136),  4*  on  this  boundary  was  taken  as 
4/  =  0  . 

Both  Eqs.  A-8  and  A-9  are  satisfied  by  the  form,  Eq.  A-10  .  When 
combined  with  Eqs.  A-3  and  A-4  to  eliminate  the  fictitious  values  (|>4  and 
4^4,  and  simplified  by  use  of  Eq.  A-10,  Eq.  A-7  becomes 

h 


4>2+4'2“  \  l+[l"'’^] 


a 


><j>o  +  horllz 


1  + 


:  0 


Eq,  A-10  and  A- 11  were  used  for  the  determination  of  stress 
functions  on  this  boundary. 

Similar  manipulation  of  equations  is  necessary  On  all  boundaries  to 
achieve  the  most  useful  form  of  the  boundary  conditions. 

Boundary  q  =  H 

Eq,  23  when  combined  with  Eq.  17  can  be  written  as 
:  p(nj  +  HZI  . 

which  can  be  integrated  along  the  bovmdary  to  give 


4i  s 


£ 


[113+ HZ] 


(A-5) 

(A-6) 


(A-7) 

(A-8) 

{A.9) 


(A-10) 


(A-11) 


(A-12) 


(A.  13) 
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Eq.  E4  when  combined  with  Eq,  18  and  put  into  finite  difference  form 
yields 


-  ‘Iji  ■  0  * 

which  when  combined  with  Eqs.  A-3  and  A-4  to  eliminate  the  fictitious 
functions  <t>i  and  i|ji  gives 

■2~  +  <j>2  +  <|>4  -  2(j)o  +  2^^3  T  Ekjjio  -  0  . 


On  the  boundary  q  =  H,  the  stress  functions  were  determined  with 
the  use  of  Eqs.  A- 13  and  A- 15 
Boundary  pel 

Eqs.  25  and  9  give  one  stress  function  directly  in  terms  of  the 
other, 


= 


1 

1  -V 


d;  +  _i-  .iinL 

^  1  +  V  1  +  V 


The  boundary  condition  given  by  Eq.  26, 


u  =  0 

n 


can  be  expressed  as 


but  with  neglect  of  rigid  body  displacements.  However  proper  applica¬ 
tion  of  Eq.  10  for  calculation  of  u*  will  alleviate  this  difficulty. 

Hence  the  boundary  condition  Eq.  26  can  be  expressed  equivalently 
as,  with  the  use  of  Eq.  14  , 


qZ  =  0 


or,  in  finite  difference  form, 


I 


4^2  -  4'4  +  V  [c|)2  *  <|)4]  -  y--—  [6  -f-  11Z]  2  0  . 

When  the  fictitious  functions  1^2  and  (1>2  are  eliminated  from  Eq.  A-20  by 
use  of  Eq.  A-3  and  A-4,  there  results 

+  ']r'3  +  4'4- 


3  -4v 

TTTT 


4'0  +  V<|)4  +  I  ^  [  2  -  h] 


Eqs.  A- 16  and  A-21  were  used  for  the  determination  of  <|>  and  ij;  on 

this  boiuidary,  * 


(A- 14) 

{A-15) 


(A- 16) 
(A-17) 

(A.18) 

(A-19) 

(A-20) 

f 

(A.21) 
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Boundary  r|-  Q 
1.  Symmetrical  Problems 

From  Eqs.  27  and  28,  the  symmetry  conditions  can  be  reduced 


to 

4»3  ■  4>i 

(A.22) 

and 

4)3  s  4)1  , 

(A-23) 

which  give,  by  elimination  of  (|>3  and  4a  the  fictitious  values  Eqs.  A-3  and 
A-4,  the  equations  which  apply  to  the  stress  function  on  this  boundary. 


[4*4  -  4>2l  +  24>i  +  <1>2  +  414  -  44>o  a  0 
and 


2^  [  4^4  -  'I'zl  +  24ii  +  4<2  +  4*4  -  44^  -  24)1  +  :  0  • 

2,  Unsymmetrical  Problems 

(a)  When  this  end  is  bonded  to  a  rigid  case,  Eqs.  29  and  30 
apply.  Eq.  29.  with  Eq.  9,  gives  a  direct  relation  between  the  stress 
functions 

2 

4/  +  [  1  -V  1<^  +  z  0 


which  simplifies  for  ti  :  0  to 

-2  e 

4j  +  [  1  -  v]4>  +  -jtj-jj-  :  0  .  ^ 

Eq.  30  is  identically  satisfied  by  the  definition  of  Eq.  10,  if 
f(  p)  ;  0  and  thus  yields  no  additional  boundary  relations.  However,  if 
u^  5  constant  when  q  r  constant,  it  is  implied  that 

"5^  =  °  • 

This  relation  is  not  useful  applied  to  Eq.  30.  However  from  the  stress- 
strain  relations 


2  [  1  +  v] 


a 

pq 


and  by  definition 

8u  8u 

C  -  p  4.  .......  -tl 

pq  oq  op 


2[  1  -f  v]  94) 
p  5q 


( A- 24) 
(A-25) 


(A.26) 

(A-27) 


(A-28) 


(A-29) 


-19- 


Combining  the  last  three  equations,  there  results 

2[l  +  v]  M 

H  "  P  M 
du 

or,  on  expressing  terms  of  the  stress  functions, 


2[  1  +  v]  =  [  I  +  v]  +  U  “  v] 


-  vp®Z 


ir\  -  ^  -  |^8fj  ‘  ‘on 

which  can  be  written  in  finite  difference  form  as 

i{ii  -  -  [  1  -  v]  t4>i  Z  :  0  . 

Elimination  of^/j  and<t»3,  at  the  fictitious  node  from  Eq.  A-32  using 
Eqs.  A-3  and  A-4  gives 

<l»l  -  2(|>o  -  24(1  +  [  1  -  2v]  ^(|)l  +  (t>2  +  <t>4  ■“  [<l»4  -<|>2l 

-  [6  +  vhZ]  s  0  . 

Thus  Eqs.  A-27  and  A-33  are  the  most  useful  forms  for  the  deter¬ 
mination  of  stress  function  on  this  boundary  for  this  end  condition. 

(b)  For  a  free  end  but  with  unsymmetrical  loading,  the 
boundary  conditions  are  the  same  as  for  the  end  q  =  H  and  may  be  reduced 
to  the  form 


-  a' 
- 


or,  since  q  =  0 


4(  = 


L 


T 


[Oi+qZ] 


n 


and 


^  [4>4  “  4*2]  4*2  +  4*4  ~  24>o  +  24(1  -  24(0  ~  0 


Eqs.  A-35  and  A-36  suffice  for  the  boundary  equations  on  this 
boundary  for  this  end  condition. 

Corner  Nodes 

At  the  corner  nodes  special  consideration  must  be  made  since 
boundary  conditions  for  both  boundaries  must  apply.  Since  this  consists 
mainly  of  algebraic  manipulation  of  equations  previously  derived,  details 
will  be  omitted  and  only  the  essential  steps  will  be  mentioned. 


{A-30) 


(A.31) 


(A.32) 


(A-33) 


(A- 34) 

(A-35) 

(A-36) 
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Corner  ti  ;  H,  p  a  a 

Boundary  conditions  Eqs.  22,  23,  and  24  can  be  satisfied  by  direct 
calculation  ofiji  using  Eqs.  A- 10  and  A- 13.  Eq,  21  will  be  satisfied  by 
applying  Eq.  A- 1 1  directly. 

Comer  p  s  1,  t|  s  H 

The  boundary  conditions  at  this  node  are 

u  :  0  ^0  =  0 

P  ® 

or,  equivalently, 

u_^  t  0  €  5  0 

*1  n 


(T  z  ■”  Us 

*1 


Applying  stress-strain  relations  of  the  form 


the  state  of  stress  and  strain  at  this  corner  is  completely  determined 
to  be 


Ha 

»l 

otq  =  -  ns 


(A-37) 


(A-38) 

(A.39) 


a 

p 


€ 

P 


and 


,L:i  ^ 

V 

1 1  +  v]  -  [  1  -  zv]  n,} 


(A-40) 


Unfortunately,  it  is  not  possible  to  satisfy  all  of  the  relevant  finite 
difference  equations  at  this  node  in  a  normal  manner.  Difficulty  with 
convergence  of  analytical  solutions  in  similar  instances  has  been  recorded 
(  see,  for  instance,  Pickett  [7]).  The  practical  reason  for  this  difficulty 
is  that  this  is  a  region  of  rapidly  varying  stresses*.  Finite  difference 
approximations  are  valid  only  when  functions  vary  slowly  with  respect  to 
the  interval  length.  It  would  be  expected  that  the  finite  difference  solution 


*The  possibility  that  a  singularity  exists  at  this  point  has  postulated 
by  many  but,  to  the  writer' s  knowledge,  has  not  been  demonstrated. 
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would  converge  to  the  correct  solution  as  the  mesh  size  approaches  zero. 
However,  for  any  finite  grid  size,  the  stresses  in  the  vicinity  of  this 
corner  cannot  be  accurately  determined.  The  effect  of  this  error  should 
be  negligible  several  nodes  from  this  corner,  this  was  confirmed  by 
obtaining  the  solution  for  one  particular  problem  with  several  mesh  sizes, 

For  the  actual  determination  of  <|>  and  at  this  corner,  Eqs,  A- 13 
and  A- 16  are  used.  This  allows  the  determination  of  and  i^i  directly. 
The  fictitious  nodes  (to  be  used  for  calculation  of  stresses,  strains  and 
displacements  as  explained  below)  are  determined  from  the  other  two 
boundary  conditions  and  one  of  the  two  field  equations.  The  other  field 
equation  is,  in  general,  unsatisfied. 

Corner  n  c  0»  p  s  1 

1.  Symmetrical  Problems 

For  the  symmetry  conditions,  Eqs.  A- 16  and  A-21  are  used 
'with  the  substitution 

>  4<i  • 

2.  Unsymmetrical  Problems 

(a)  For  the  bonded  end  at  this  node,  the  conditions 

u„:u  :€-:£=€:€  cO 

n  p  0  p  q  pr 

apply.  These  conditions  are  satisfied  by  using  Eq.  A-16  and  A-21  with 
v{»3  ■  in  Eq.  A-21. 

{ b)  For  a  free  end  with  unsymmetrical  loading,  this  node  is 
treated  the  same  as  the  point  p  r  1,  q  :  H  . 

Corner  g  5  0,  p  zta 

1.  Symmetrical  Problems 

For  the  symmetry  condition  Eqs.  A- 10  and  A- 11  are  used 
unchanged. 

2,  Unsymmetrical  Problems 

(a)  For  the  bonded  end  Eq.  A- 10  and  Eq.  A-27  are  used  with 
the  same  type  of  error  as  noted  for  the  node  p  a  1,  q  =  H. 

(b)  Eqs.  A- 10  and  A- 11  also  apply  for  unsymmetrical  loading. 
Method  of  Solution 

If  the  finite  difference  equations  are  written  at  each  node  of  an  MxN 
mesh,  the  stress  functions  at  each  node  may  be  found  by  simultaneous 
solution  of  the  resulting  set  of  {  M  +  I)  (  N  +  1)  equations.  In  practice  this 
may  involve  sets  containing  on  the  order  of  1000  equations.  A  number  of 


(A-41) 


(A-42) 


III 


methods  may  be  chosen  to  achieve  a  solution,  but  the  common  method 
used  in  elasticity  problems  is  an  iterative  procedure,  jbpulariaed  by 
Southwell  [  I3j,  known  as  the  relaxation  method. 

The  Seidel  iteration  process  using  systematic  over- relaxation  was 
used  [  14]  to  obtain  the  values  of  the  stress  functions.  New  values  of  a 
stress  function  <|).  at  nodal  point  i  for  the  k+  1  iteration  are  found  from 


<|>/  =  (J.  +  ^R.  , 


( k) 

where  4..  is  the  stress  function  value  of  the  i  node  for  the  kth  iteration. 
0  18  the  over-relaxation  factor,  and  R.  is  the  residual  of  the  finite  differ¬ 
ence  equation.  The  speed  of  convergence  is  a  strong  function  of  ^  which 
must  be  in  the  range  1  <  ^  <  2.  Values  of  I.  5  to  1. 7  for  /3  were  found  to 
be  near  optimum,  but  the  optimum  value  varies  with  the  input  parameters 
of  each  problem.  Starting  with  an  estimate  of  zero  for  each  stress 
function  value  for  a  4x32  mesh  with  a  s  1/3,  some  200  to  300  iterations 
on  each  stress  function  were  required  with  jS  r  1.  7  to  obtain  the  solution 
for  the  stress  functions  with  the  residuals  reduced  to  0.000001  of  the 
largest  stress  function  value. 

Because  the  stress  function  equations  are  coupled  at  the  boundaries 
and  the  interior,  instability  problems  can  arise  under  certain  conditions. 

For  example,  the  following  process  was  found  to  work  satisfactorily  for 
symmetrical  loading: 

1.  Initially  assume  all  stress  function  values  to  be  zero. 

Make  one  iteration  correcting  the  values  of  ^  on  the  boundaries 
p  :  tt  and  fj  :  H  with  fi  =  1.  Corner  nodes  may  be  treated  as  part  of 
either  of  the  intersecting  boundaries. 

3.  Relax  the  interior  values  and  boundary  values  at  ?i  =  0  of  <}» 
usingg:  1.7.  iterating  until  Eqs.  A.3  and  A.24  are  satisfied,  i.e..the 
absolute  value  of  the  largest  residual  is  less  than  some  specified  value. 

^  4.  Make  one  iteration  correcting  the  values  of  +  on  the  boundary 
p  :  1.  4^  is  known  exactly  on  the  boundaries  ft  =  o  and  q  s  H. 

5.  Relax  4-  at  q  =  0  and  in  the  interior  until  Eqs.  A-4  and  A.25  are 
satisfied  to  some  specified  degree  as  in  step  (  3)  above. 

6.  Calculate  4>  on  boundary  p  =  1  by  Eq.  A- 16  and  return  to  8tep(2). 
The  problem  was  considered  relaxed  to  a  specified  degree  when  the 

absolute  value  of  all  residuals  was  less  than  a  specified  number  d.  and 
the  sum  of  all  residuals  was  less  than  a  specified  number  D,  i.e., 


(A-43) 


(A-44) 


]CR-|  <  D  . 

The  useful  value  of  O  is  relatedto  both  d  and  the  minaber  of  node 
points.  0  was  estimated  by  the  equation 
r.  (mnI  . 


(A.45) 


where  M  is  the  number  of  increments  in  the  axial  direction  and  N  is  the 
number  of  increments  in  the  radial  direction. 

The  relaxation  process  was  found  to  go  faster  if  the  value  of  d  was 
initially  chosen  large  (say  0,001  of  the  largest  stress  function  value)  until 
the  problem  was  relaxed  completely;  then  d  was  reduced  by  an  order  of 
magnitude,  and  the  problem  further  relaxed  to  completion  with  respect  to 
the  new  value  of  d;  then  the  process  was  repeated  until  the  desired  degree 
of  relaxation  was  achieved. 

Calculation  of  Stresses,  Strains  and  Displacements  From  the  Stress 
Functions 

The  stresses,  strains,  and  displacements  are  calculated  from  the 
stress  functions  by  writing  the  finite  difference  forms  of  Eqs,  9  through 
18.  Since  these  equations  involve  no  derivatives  higher  than  the  first 
this  process  is,  in  general,  straightforward  and  results  in  high  accuracy. 
In  practice,  two  factors  complicate  the  calculation  of  these  quantities. 
First,  the  derivatives  normal  to  a  boundary,  when  written  in  terms  of  the 
central  finite  difference  formulae,  involve  fictitious  nodes;  or,  when 
written  in  terms  of  forward  or  backward  differences,  involve  nodes  lying 
only  on  one  side  of  the  boundary.  The  latter  method  was  found  to  introduce 
unduly  large  errors,  and  recourse  was  made  to  the  calculation  of  the 
fictitious  nodes  needed  to  express  the  normal  boundary  gradients  in  terms 
of  central  difference  formulae.  Secondly,  the  axial  disj^acement,  llq.  10, 
involves  integrals  in  the  axial  or  q  direction  of  derivatives  with  respect  to 
the  radius  p.  It  was  found  that  the  application  of  Simpson' s  rule  for  this 
integration  was  satisfactory  at  even  numbered  nodes,  although  it  was  quite 
complicated  to  program  for  a  computer.  The  accuracy  of  the  axial  dis¬ 
placement  at  odd  nodes,  where  trapezoidal  integration  was  necessary  for 
one  increment,  was  not  always  satisfactory. 
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Computation  Facilities 

This  problem  has  been  programmed  for  the  IBM  7090  computer  at 
the  NASA  Marshall  Space  Flight  Center.  Typical  machine  times  for  pro¬ 
blems  of  this  type  vary  from  about  one  minute  for  problems  with  50  nodes 
to  20  minutes  for  500  nodes.  These  times  are  also  greatly  dependent  on 
Poisson' s  ratio,  v,  and  the  radius  ration  a. 
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Fig.  0  Efftd  of  Poisson’s  ratio,  w,  on  profile,  propellant  bonded  at  outer 
periphery,  both  ends  free,  under  uniform  shrinkoge 
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Fig.  9  Effect  of  radius  ratio,  ic  ,  on  profile,  propellant  bonded  at  outer 
periphery,  both  ends  free,  under  uniform  shrinkage  X*l  ,  ^'•O.S 
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Radiol  displacement  at  Inner  radius  in  middle  of  cylinder  of  finite  length  as  p>er- 
centoge  of  rodiol  displacement  in  equivalent  cylinder  of  infinite  length,  propellant 
bonded  at  outer  periphery  ,  both  ends  free ,  under  uniform  shrinkoge  vo.49 
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Length-to-diameter  ratio,  \ 


Radial  displacement  at  inner  radius  in  middle  of  cylinder  of  finite 
length  os  percentage  of  rodiol  displacement  in  equivolent  cylinder 
of  Infinite  length,  propellant  bonded  ot  outer  periphery,  both 
ends  free,  under  uniform  shrinkage  t'»o.45 
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Rg.  16  Radiol  displocemenf  at  Inner  radius  In  middle  of  cylinder  of  finite  length  os 
percentage  of  radial  displacement  in  equivalent  cylinder  of  infinite  lenath 
propellant  bonded  ot  outer  periphery,  both  ends  free  ,  under  liS 
shrinkage  v'O.so 
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Dimensionless  circumferential  stress  at  inner  rcxlius  in  middle  of  cylinder ,  propellanf  bcwided  at  outer 
iphery  ,  both  ends  free  ,  under  uniform  shrinkage  *'»o.5 
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Length-to-diameter  ratio,  A 


Rg  19  Dimensionless  clrcumfcrenild  stress  at  Inner  radius  in  middle  of  cylinder, 
propellant  bonded  at  outer  periphery,  both  ends  free  ,  under  uniform 
shrinkage  «/*o.45 


Length-to-diameter  ratio,  \ 


g.  dV  D  mensionicss  qaIqI  stress  at  inner  radius  in  middle  of  cylinder,  propell 
bOTded  at  outer  periphery,  both  ends  free,  under  uniform 
shrinkage  i/-o.45 
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Rg  21 


DitDenslonless  circumferential  stress  at  inner  radius  In  middle  of  cylindei 
propellant  bonded  ot  outer  periphery,  both  ends  free,  under  uniform 
shrinkage  v'o.as 


Length-to-diameter  ratio,  \ 


Dimensionless  axial  stress 


Fig-  22  Dimensionless  qxIqI  stress  at  Inner  rodius  in  middle  of  cylinder, 

propellant  bonded  at  outer  periphery,  both  ends  free,  under  uniform 
shrinkoge  v-o.35 
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Rg.  23  Deformation  of  propellant  grain,  bonded  at  outer  periphery,  both  ends  free,  under  oj 
»f  ‘3,33  ,  X«  2  ,  v»0.5 
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"ig,  24  Deformation  of  propeltont  grain,  bonded  at  outer  periphery  and  left  end,  under  qxIoI 
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ends  free ,  under  qxIqI  occelerafion 


25  Dlmmsiontess  3lr«5  corttourj  5*6  In  propellant  groin,  bonded  at  outer  periphery  ond  left 

under  oxld  occelerdlon  k.3  ,  X»2.67  ,  i/.  0,5 
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DISPLACEMENT  ANALYSIS  BY  ENERGY  METHODS  FOR 
CURE  SHRINKAGE  OF  CASE-BONDEt)  SOLID  PROPELLANT  GRAINS 

OF  FINITE  LENGTH 


Gerald  F.  Gillis 

The  use  of  energy  methods  for  the  approximate  solution  of  elasticity 
problems  is  well  established  for  certain  types  of  bodies  for  which  at 
least  one  overall  dimension  is  small  in  comparison  with  the  others,  e.g., 
beams  and  plates.  The  investigation  of  the  stress  state  and/or  displace¬ 
ment  state  in  bodies  in  which  all  overall  dimensions  are  of  the  same 
magnitude  has  received  little  attention  by  energy  methods,  evidently 
because  the  simplicity  evident  in  plate  and  beam  problems  is  no  longer 
present.  This  report  has  served  as  an  exploratory  study  of  the  solution 
of  problems  of  the  latter  type.  The  particular  problem  chosen  for  the 
investigation  was  selected  because  of  the  availability  of  both  numerical 
analysis  and  experimental  results  with  which  to  verify  the  solution. 

PROBLEM  DESCRIPTION 

Strains  are  induced  in  polymeric  materials  as  a  result  of  thermal 
expansion  and  cure  shrinkage.  Due  to  the  similarity  of  the  two  effects, 
they  may  be  combined  into  a  dimensional  change  function,  A,  defined  as 

A  =  6  +  nT  .  (  1) 

where 

6  =  linear  cure  expansion  {positive  for  expansion  and  negative  for 
shrinkage)  , 

O'  =  linear  coefficient  of  thermal  expansion  of  the  cured  material,  and 

T  =  temperature  relative  to  some  reference. 

With  the  notation  of  Eq.  1,  (  €  -  A)  represents  strain  due  to  stress 
whereas  €  represents  overall  dimensional  change.  When  there  are  no 
applied  loads  or  restraints,  €  r  A  and  the  stress  induced  strain  is  zero. 

If  the  effects  of  the  thermal  expansion  and  cure  shrinkage  are  such  that 
A  is  zero,  then  the  stress  induced  strain  is  €,  as  it  is  usually  defined. 

Stresses  are  induced  in  case-bonded  propellant  grains  by  cure 
shrinkage  and  by  uniform  temperature  changes,  the  latter  being  the 
result  of  the  difference  in  the  coefficients  of  thermal  expansion  of  the 
case  and  grain.  In  determining  the  stresses  caused  by  the  above  effects, 
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T  would  normally  be  the  temperature  relative  to  the  temperature  at 
which  no  stresses  are  induced  by  differential  thermal  expansion  of  the 
case  and  grain. 


The  purpose  of  this  analysis  was  to  determine  by  energy  methods 
the  displacements  of  the  propellant  grain  associated  with  cure  shrinkage 
only.  However,  the  results  are  general  enough  to  include  thermal  expan¬ 
sion  resulting  from  uniform  temperature  changes  if  the  coefficient  of 
thermal  expansion  of  the  case  is  negligible  compared  to  that  of  the  grain. 

The  model  which  was  analyzed  consisted  of  an  elastic,  homogeneous, 
isotropic,  incompressible  propellant  grain  in  the  form  of  a  hollow,  right, 
circular  cylinder  bonded  on  the  outer  surface  to  a  rigid  case  wall.  The 
applicability  of  the  classical  infinitesimal  theory  of  elasticity  was 

assumed. 

Cylindrical  coordinates  were  used  with  the  origin  located  at  the 

centroid  of  the  cylinder 

(  Fig.  1) .  The  analysis  was  made  only  for  the 

half  of  the  cylinder  extending  from  z  :  0  to  z  :  i  since  it  is  symmetrical 

about  the  origin. 

The  boundary  conditions  for  the  described  model  are 

a) 

u  =  0,  r  =  b  , 

b) 

w  =  0,  r  =  b  , 

c) 

=  0,  r  =  a  , 

d) 

=  r  =  a, 

e) 

f) 

^2  =  0*  z  =  i  , 

g) 

^rz  ~  z  =  0  ,  and 

h) 

w  =  0,  Z  r  0  , 

where 

u 

=  radial  displacement  , 

w 

=  axial  displacement  , 

T 

=  shear  stress, 

a 

=  normal  stress  , 

a 

=  inner  propellant  radius  , 

b  . 

=  outer  propellant  radius, 

1 

=  one-half  the  length  of  the  propellant  grain  , 

z 

=  axial  direction  , 

r 

=  radial  direction  ,  and 

e 

=  circumferential  direction  . 

) 
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Boundary  conditions  (a)  and  (b)  are  a  consequence  of  the  rigid  case 
assumption.  The  restriction  to  loads  caused  by  cure  shrinkage  requires 
( c)  I  ( d)  I  { e) ,  and  ( f) .  Conditions  ( g)  and  ( h)  result  from  symmetry 
considerations.  The  axial  symmetry  of  geometry  and  loading  require 
that.  =  r^^*Tg^:0,  where  v  is  circumferential  displacement. 

ANALYSIS 

An  approximate  displacement  solution  for  the  described  model  was 
obtained  by  applying  the  Rayleigh-Ritz  method  to  the  principle  of  mini¬ 
mum  potential  energy*.  This  was  accomplished  by  assuming  the  dis¬ 
placements, in  equation  form  with  undetermined  constants  such  that  all 
displacement  boundary  conditions  were  satisfied  and  the  displacements 
were  continuous.  These  displacements  were  substituted  into  the  potential 
energy  equations  and  the  undetermined  constants  were  evaluated  by  mini¬ 
mizing  the  potential  energy.  The  approximation  was  improved  by  increas¬ 
ing  the  number  of  undetermined  constants  and  repeating  the  minimization 
procedure  until  the  potential  energy  approached  a  stationary  value. 
Potential  Energy 

The  boundary  conditions  for  the  model  state  that  the  surface  tractions 
are  zero  at  r  s  a  and  z  s  i  i  and  that  us  vswsOatrsb.  Consequently, 
for  negligible  body  forces,  the  potential  energy  as  defined  by  Eq.  A-8  is 
equal  to  the  strain  energy.  That  is, 


where 

^  =  volume  of  the  body  , 

E  =  potential  energy  , 

Jr 

V  =  strain  energy  in  the  body,  and 
=  strain  energy  per  unit  volume  . 

The  strain  energy  per  unit  volume  for  a  body  subjected  to  a  linear  cure 
expansion  6  is  defined  by  Eq.  A- 5  as 


2V 


o  “  ^]  +  O'  -  6]  +  T  y  +  T  „  Y  ,  +  T  V 
6  9^  Z*-  z  J  rz  ’rz  rS^re  zQ'ze 


(3) 


*The  Appendix  contains  a  discussion  of  energy  methods  which  includes 
the  equations  and  limitations  of  the  principle  of  minimum  complementary 
energy  as  well  as  the  principle  of  minimum  potential  energy  and  the  pro¬ 
cedure  for  applying  the  Rayleigh- Ritz  method  to  either. 
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UBing  Hooke' s  law  in  the  (orm 

OTq  -  J'^"y  ®  **■  2G[Cq-6] 

-  tit  0 

where 

G  c  shear  modulus  , 

0  s  CTj.  +  ffg  +  ,  and 

V  s  Poisson'  a  ratio  , 

Eq.  3  can  be  written  as 

V 


and 


(4) 

(5) 

(6) 


2V 


o  =  e  [  e  -  36]  +  2G[€  /  +  -  26e  +  36^] 


+  G[Vj./+Vi,e^+\e.^l  * 

where  e  s  +  €q  +  .  For  an  incompressible  material,  e  -  36  s  0  . 

Also,  as  stated  above,  T^g  =  Tq^  z  0  due  to  symmetry.  Hence  Eq.  7  can 
be  written  as 

V  :  G[€  *  +  e.H  €  ^  -  36*]  +  £  V  2  , 
o  '•  r  0  z  ^  I  'tz  ' 

Approximating  Functions  for  the  Displacements 

Parr  [  6]  *  has  shown  that  the  radial  displacement  in  a  propellant 

grain  in  the  form  of  a  hollow  circular  cylinder  of  infinite  length,  bonded 

to  a  rigid  case  and  subjected  to  cure  shrinkage,  is 

u2 


[  1  +  v]  a*  r  6 

US-  -V - i - 

a*  +  [  1  -  2v]b* 


-  1 


where  -  6  is  the  linear  cure  shrinkage.  For  an  incompressible  material 
V  :  0.  5,  and  Eq.  9  simplifies  to 

^  =2  6[p  -  p"^] 

where  p  :  r/b  . 

The  radial  displacement  for  a  cylinder  of  finite  length  can  be  assumed 
to  consist  of  the  infinite  cylinder  solution  plus  a  perturbation  term, 

E  :!  8[p-p'‘)  +5fr(p)  G(c)  , 

where  ^  r  j  and  F(  p)  and  G(^)  are  as  yet  unspecified  functions. 


(7) 


(8) 


(9) 


(10) 


(11) 


‘Numbers  in  brackets  designate  references  at  the  end  of  the  report. 
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Forming  an  assumed  expression  for  the  axial  displacement  in  a 
similar  manner  results  in  only  the  perturbation  term  since  w  s  0  in  the 
infinite  cylinder  solution.  Hence, 

7-6f{p)  g(&)  . 

where  f(  p)  and  g(  jj)  are  again  as  yet  unspecified  functions. 

The  restriction  to  an  incompressible  material  requires  that 


+  €  +  €  = 
6  z  r 


u  aw  8u 

7  ^ 


.  1  av  -- 

7  75  =  “ 


Substitution  of  Eqs.  11  and  12  into  Eq.  13  with  v  a  0,  yields  the 
conditions 

g'(;)  :\G(t) 


where  A  is  an  undetermined  constant  which  can  be  evaluated  by  minimizing 
the  energy.  In  this  analysis,  \  will  be  taken  equal  to  unity  for  simplicity. 

Experimental  measurements  of  axial  and  radial  displacements  made 
by  the  Applied  Physics  Group  of  this  Division  on  a  6C2-33  STM*  (  Figs. 

2  and  4)  suggest  using  a  parabola  for  f(p)  of  the  form 

f(p)=4[l-p^]  , 

where  the  factor  of  4  was  introduced  to  obtain  a  more  convenient 
expression  for  F(p)  .  Substitution  of  Eq.  16  into  Eq.  15,  integration  of 
the  resulting  expression,  and  application  of  boundary  condition  (a)  yields 

F( p)  =  [  p  *  +  p^  -  2p] 

These  radial  and  axial  displacement  measurements  further  suggest 
the  possibility  of  employing  polynomials  for  G(  and  g{  .  Furthermore, 
symmetry  of  the  radial  displacements  about  the  9  -  z  plane  requires  that 
a  polynomial  representation  of  G(  ;)  contain  only  even  powers  of  Con¬ 
sequently,  it  was  assumed  that 


2i-  1 


*602-33  STM  designates  a  static  test  motor  (STM,)  having  a  case- 
bonded  propellant  grain  in  the  form  of  a  hollow,  right,  circular  (  C) 
cylinder  33  in.  long  with  an  O.D.  of  6  in.  and  an  l.D.  of  2  in. 
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,21-.  2 


G(t)  =/(t)  [2i-l]  C2..1  C 

Substitution  of  Eqs.  16, 17, 18,  and  19  into  Eqs.  11  and  12  yields 


(19) 


^  s  I  6  [  P  -  p’^]  +  6  [ p‘*  +  p’  -  2p]  )  [  2i  -  1]  C2J  _  j  & 


2i-2 


fa.li I  III 

i:l 


(20) 


and 


w 

T 


46  [1 


-4 


'2i-  1 


2i  -  1 


is  1 

Eqs,  20  and  21  satisfy  all  displacement  boundary  conditions  and  are 
continuous . 


An  e>qponential  axial  variation  was  also  considered  in  an  effort  to 
determine  which  type  of  function  gives  the  best  approximation  with  the 
least  expenditure  of  time.  The  solution  yielding  the  smaller  potential 
energy  should  be  the  better  approximation. 

The  displacements  with  the  exponential  axial  variation  were  assumed 
to  be 


and 


w_6B6bri  )  kjfr.  il 

j  =  [  1  -  P  1  exp  -^  [  1  -  ^]  j 

where  B  and  k  are  the  undetermined  constants  and  k  is  positive  and  real. 
These  displacements  are  continuous  and  satisfy  all  displacement  boundary 
conditions  with  the  exception  that  the  assumed  form  for  w  does  not  exactly 
satisfy  the  boundary  condition  that  w  ■  0  at  z  s  0.  This  boundary  condition 
will  be  approximately  satisfied  if  the  exponential  term  at  z  =  0  is  small 
relative  to  its  value  at  z  s  i.  Eqs.  22  and  23  were  written  to  apply  only 
for  the  region  0  :£  z  :<  1  with  the  knowledge  that,  if  a  good  approximation 
could  be  obtained  for  this  half  of  the  motor,  symmetry  about  the  origin 
will  require  the  other  half  to  exhibit  the  same  behavior. 

Substitution  of  the  displacements  given  by  Eqs.  20  and  21  and  Eqs. 

22  and  23  into  Eq.  8  and  the  resulting  expressions  into  Eq.  2  yields 
potential  energies  for  the  half-length  of 


(21) 


(22) 

(23) 


6*jrG|b^ 


i  ■  1  j  s  1 
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[2Y.+  32V.1 


[2i+2j-5]2‘ 


[2i-lll2j-2]  +r2i.l][2i-2] 


b*[2i+  2j  -  1] 


Ji  +  2j  - 


+  3  (3v2  -  Zys)  +  12^3 


[16Y4] 


t  +  im.  +  +  9  L  .2  _] 

b^nGlh^  2i  1^  k  4  j^3  J  ki  3  J 

for  the  polynomial  and  exponential  solutions,  respectively,  where 


Yi  =  4ri^  -  I  ti6  _  ir^  +  ^^-2  _  ^  , 

V!  — id'  -  ’ll  . 

VJ  =2l  1  -1*+  ’1*-1-1  , 

Y4  =  2  -j-n* +  >1* , 


Y5  =  j  [1  -’1'^  , 

Y6  =  -  -  In  9  -  2  +  2V  -  ^  9*  +  ^  Il‘  , 

Yt  =  7  1 -  3n‘  -  2  -  ‘*’1*  ' 

Y8  =  [I  - ’t'*]  »  and 
Ti=  a/b  . 

A  sufficient  number  of  equations  for  determination  of  the  constants  may  be 
obtained  by  taking  the  derivative  of  the  above  potential  energies  with 
respect  to  each  constant  and  equating  the  derivative  to  zero.  The  process 
is  much  too  cumbersome  to  be  carried  out  here. 


. . . mill- 
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NUMERICAL  RESULTS 

Numerical  calculations  were  made  for  a  propellant  grain  with  inner 
radius  a  *  1  in.,  outer  radius  b  =  3  in.,  half  length  i  =  16  in.,  and  a 
linear  cure  shrinkage  -  5  =  0. 00428  in.  /in. 

Minimization  of  the  potential  energy  as  defined  by  Eq.  24  retaining 
only  three  constants  (three  terms)  in  the  assumed  polynomial  form  for 
the  displacements  yielded 
Cl  =  0.047897  , 

C3  =  -  0.068692  , 

C5  =  0. 289720  ,  and 
E 

— 2 —  =  30.9525  . 

ff6*Gib* 

The  following  values  were  obtained  when  the  number  of  constants  was 
increased  to  four. 

Cl  =  0.005841  , 

C3  =  0.239486  , 
cs  = -0.2773.36  , 

C7  =  0, 303037  ,  and 
E 

- -  30.8979  . 

7r6^Gib* 

The  addition  of  the  fourth  constant  resulted  in  an  energy  decrease  of  only 
0,  18% .  Consequently,  additional  constants  would  not  be  expected  to 
significantly  improve  the  approximation. 

Minimization  of  the  potential  energy  as  defined  by  Eq,  25  for  the 
exponential  form  of  the  displacements  yielded 
k  =  0.948515  , 

B  =  0.911945  ,  and 
E 

-  =  30.  8979  . 

r6*Gib^ 

With  the  above  constants,  the  value  of  the  exponential  term  in  Eq.  23  at 
z  =  0  is  two  orders  of  magnitude  smaller  than  its  value  at  z  =  i;  hence, 
the  approximation  tow  =  0atz  =  0i8  reasonable. 

The  radial  displacement  along  the  inner  surface  of  the  propellant  for 
each  of  the  above  solutions,  for  a  solution  obtained  numerically  using 
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relaxation  techniques^  and  from  measured  values  is  shovm  in  Figs.  Z  and 
3.  The  associated  axial  displacements  at  z  =  |  are  shown  in  Fig.  4. 

It  should  be  noted  that  the  numerical  values  for  the  above  constants 
are  independent  of  the  cure  shrinkage  but  are  functions  of  the  ratios  ^ 
and  j  ;  hence,  they  are  valid  for  any  motor  where  ^  and  j  =  ^  . 


DISCUSSION 

As  was  previously  mentioned,  additional  undetermined  constants  are 
added  to  the  approximating  functions  until  the  energy  approaches  a  sta¬ 
tionary  value.  The  fact  that  a  stationary  value  is  obtained  does  not 
necessarily  imply  a  close  approximation  to  the  exact  solution,  A  sta¬ 
tionary  value  merely  assures  that  the  solution  is  about  as  good  as  can  be 
obtained  from  the  chosen  function  forms.  If  the  approximating  functions 
are  poor,  the  approximation  may  be  incapable  of  approaching  the  true 
minimum  energy  or  the  convergence  to  the  true  minimum  may  require  a 
large  number  of  constants.  Experience  in  applying  energy  methods  and 
experimental  data  are  valuable  assets  in  choosing  functions  with  the  best 
chance  of  a  close  approximation. 

The  solutions  obtained  above  could  be  improved  by  an  iteration  pro¬ 
cedure.  For  example,  definite  expressions  were  assumed  for  f(  p)  and 
F(  p)  and  the  undetermined  constants  which  were  included  in  G(  and 
g(  £,)  were  evaluated  by  minimizing  the  potential  energy.  With  these 
expressions  for  G(  and  g{  ,  undetermined  constants  could  then  be 
introduced  into  F{p)  and  f(  p)  and  evaluated  by  energy  minimization.  The 
process  is  then  repeated  as  many  times  as  desired.  For  the  same  degree 
of  accuracy,  the  iteration  procedure  has  the  advantage  over  one  which 
initially  introduces  many  unknown  constants  in  that  at  each  iteration  only 
a  few  constants  must  be  determined  through  simultaneous  solution  of 
algebraic:  equations. 

If  experimental  data  are  not  available  with  which  to  compare  results, 
upper  and  lower  bounds  on  the  energy  should  be  determined,  if  possible, 
so  that  the  maximum  deviation  from  the  true  solution  can  be  assessed. 

As  an  example,  the  problem  considered  in  this  analysis  was  concerned 

^This  solution,  obtained  by  Parr,  is  presented  in  another  section  of 
this  report. 


-56- 


with  the  cure  shrinkage  of  a  cylindrical  propellant  grain  of  finite  length. 
An  upper  bound  on  the  potential  energy  could  be  established  from  the 
infinite  cylinder  solution  (plane  strain)  obtained  by  Parr  [  6] .  A  lower 
bound  could  be  determined  by  generalizing  Timoshenko' s  polution  [  8] 
for  a  thin  circular  disk  (plane  stress)  subjected  to  a  temperature  varia¬ 
tion  symmetrical  about  the  axis  of  the  cylinder  to  include  cure  shrinkage. 

These  two  solutions  yield  upper  and  lower  bounds  of 

£ 

— £ - .=  40.0144 

7r6*Gib* 

and 

E 

— 2 -  :  4.00144  , 

jr6*Gib* 

respectively,  for  a  :  1  in. ,  i  :  16  in.,  b  =  3  in.,  and  6  :  -  0.00428 
in. /in.  A  comparison  of  these  values  with  those  given  in  the  numerical 
results  indicates  that  the  finite  cylinder  solution  is  closer  to  the  infinite 
cylinder  solution  than  to  the  one  for  a  thin  disk.  Since  the  grain  used  in 
the  example  is  relatively  long  compared  to  the  web,  it  would  be  expected 
that  the  potential  energy  would  be  near  the  upper  bound,  as  it  is.  For 
grains  which  are  short  compared  to  the  web,  it  would  be  expected  that 
the  potential  energy  would  be  near  the  lower  bound. 

CONCLUSIONS 

In  general,  the  displacements  obtained  by  energy  methods  compare 
favorably  with  measured  values  and  those  obtained  by  a  relaxation  solu¬ 
tion.  Some  departure  from  the  measured  values  was  anticipated  due  to 
the  assumption  of  incompressibility  (  v  «  0.  5) .  Present  indications  are 
that  the  propellant  from  which  the  measured  values  were  obtained  has  a 
Poisson' s  ratio  of  approximately  0.495.  However,  a  stress-strain 
analysis  is  sensitive  to  even  such  a  slight  variation  in  the  value  of  v  when 
it  is  in  the  range  0.45  to  0.50. 

The  two-constant  exponential  solution  w'as  in  Ciose  agreement  with 
the  one  obtained  using  four-constant  polynomial.  Thus,  a  better  approxi¬ 
mation  with  fewer  constants  can  be  obtained  using  the  exponential  form. 
However,  the  fact  that  the  polynomial  form  yields  linear  algebraic  equa¬ 
tions  for  the  determination  of  the  constants  is  a  factor  deserving  consider¬ 
ation.  A  solution  can  be  obtained  using  the  polynomial  form  with  more 
undetermined  constants  in  less  time  than  would  be  required  for  an 
exponential  form. 
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APPENDIX 


GENERAL  DISCUSSION  OF  ENERGY  METHODS  AND  TERMINOLOGY 


tram  Enerey  anc 


tram  Enerc 


When  a  uniform  bar  is  loaded  in  simple  tension,  the  forces  on  the 
ends  do  work  which  is  stored  in  the  body  as  strain  energy.  For  an  elastic 
element  of  unit  dimensions  shown  in  Fig.  5,  the  strain  energy  stored  per 
unit  volume  is 


a  € 
r  r 


The  strain  energy  per  unit  volume  is  equal  to  the  area  between  the 
€-axis  and  the  stress-strain  curve.  The  area  between  the  a-axis  and  the 
curve  is  called  the  complementary  strain  energy.  If  all  six  components 
of  stress  act  on  the  element,  the  principle  of  the  conservation  of  energy 
requires  that  the  work  done  be  independent  of  the  order  in  which  the 
forces  are  applied.  Otherwise  it  would  be  possible  to  load  in  one  order 
and  unload  in  another  such  that  net  work  would  be  gained  from  the  cycle. 
Consequently,  when  all  six  components  of  stress  are  considered,  the  strain 
energy  per  unit  volume  can  be  generalized  to,  in  cylindrical  coordinates, 

\  ■  ?  ["rS  +  ''e  S  +  +  "re're  +  \z'<Tz  *  "er  v]  • 

For  a  body  subjected  to  thermal  effects  and  cure  shrinkage,  the  strain 
due  to  stVess  has  previously  been  defined  as  (€  -  A),  whereas  the  strain 
due  to  dimensional  changes  is  denoted  by  €.  For  simplicity,  consider  a 
uniform  bar  subjected  to  a  uniform  temperature  change  and  then  loaded  in 
simple  tension.  The  resulting  stress-strain  diagram  is  shown  in  Fig.  6. 
The  strain  energy  per  unit  volume,  the  area  between  the  €-axis  and  the 
curve,  is 


The  complementary  strain  energy  per  unit  volume,  the  area  between  the 
ff-axis  and  the  curve,  is 

V  *  :  i  ff  ft'  f  otI  . 

0  2  r  [  r  j 

These  equations  may  also  be  generalized  to  include  all  stress  components 
and  the  cure  shrinkage.  Thus, 


(A-D 


(A-2) 


(A- 3) 


(A-4) 
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''o  =  \  |»r  +  •'e  I's-AJ  +  [e^- A]  +  ^eV,^  +  T„v„  +  (A.5) 


and 
V 


|ff,[e,  +  A]  +  ag[« j+A;]  +aje^+  Al  +  r^^yA.  (A-6) 


Principles  of  Minimum  Energy 

The  solution  of  many  elasticity  problems  can  be  obtained  employing 
energy  theorems  based,  in  part,  on  the  concepts  of  strain  energy  and 
complementary  strain  energy.  Two  such  theorems  are  the  theorem  of 


minimum  potential  energy  and  the  theorem  of  minimum  complementary 
energy.  A  detailed  discussion  of  these  theorems  is  given  by  Sokolnikoff 
[  7] .  Westergaard  [  9]  gives  an  informative  note  on  the  history  of  their 
development.  The  purpose  of  the  present  discussion  is  not  to  repeat  the 

I 

derivations  and  proofs  readily  available  in  the  literature  but  to  summarize 
the  results  as  they  apply  to  problem  solution. 

Sokolnikoff  demonstrates  that 


Ai 


where 


KmniKoii  aemonsc races  tnac  . 

/  ^  ^  Xu+  Yv+  Zw]  d|  -  ^  [Xu  +  Tv+Xw]  dS^  2:  0 


(A.7) 


=  strain  energy  per  unit  volume, 


=  summation  over  the  part  of  the  surface  on  which  the  surface 
forces  are  prescribed, 


=  summation  throughout  the  volume, 


Aj  =  an  arbitrary,  small,  continuous  change  which  is  consistent  with 
the  displacement  specified  boundary  conditions  and  associated 
with  a  variation  of  the  displacements  from  the  equilibrium  state, 

X,  y,  Z  a  the  components  of  the  body  forces  per  unit  volume  in  the  r, 

9,  and  z  directions,  respectively, 

"TV  *W  IV 

X,  Y,  Z  a  the  components  of  the  surface  forces  per  unit  area  in  the  r, 
0,  and  z  directions,  respectively,  and 

u,  V,  w  a  the  equilibrium  components  of  displacements  in  the  r,  0,and 
z  directions,  respectively. 

The  expression  in  the  parentheses  is  defined  as  the  potential  energy, 
that  is, 


f  [Xu+YvtZw]  de -J  lXu  +  Yv  +  rw)dS  . 

’  *  a 


(A-8) 
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Time,  Bq.  A-7  states  that  of  all  continuous  displacements  satisfying 

gi«n  displacement  specified  boundary  conditions,  those  which  satisfy  the 

equUabrta  conditions  make  the  potential  energy  an  absolute  minimum. 

The  above  statement  is  known  as  the  principle  of  minimum  potential 
energy. 

Sokolnikoff  also  demonstrates  that 

V  -IS  -  £  [Xu  +  Tv  +  rw]  ds]  >  0 , 

Where  ^  ^ 


®C  =  J  V  +  7v  1 7w] 


With  this  definition.  Eq.  A-9  states  that  of  all  statically  possible  combi¬ 
nations  of  stresses  that  will  satisfy  the  equilibrium  equations  and  stress 
specified  boundary  conditions,  the  one  that  corresponds  to  the  true  state 
of  stress  (the  one  that  also  satisfies  the  compatibility  equations)  wiU 
make  the  complementary  energy  an  absolute  minimum.  The  above  State- 
ment  is  known  as  the  principle  of  minimum  complementary  energy. 
Rayleigh-Ritz  Method 

An  important  application  of  the  minimum  principles  relates  to  their 
use  in  obtaining  approximate  solutions  to  problems.  In  the  Rayleigh- 
its  method,  it  is  assumed  that  a  nearly  true  solution  can  be  repre¬ 
sented  by  a  set  of  functions  involving  a  limited  number  of  constants  to  be 
determined  by  minimizing  the  energy. 

In  applying  the  Rayleigh-Rite  method  to  the  principle  of  minimum 
complementary  energy,  the  stresses  are  assumed  in  an  equation  form 

involving  undetermined  constants  such  that  the  equUibrium  equations  and 

boundary  conditions  specified  in  termi  nf  at  ^ 

..ptciuea  m  terms  of  stresses  are  satisfied.  These 

assumed  stresses  are  then  substituted  into  the  complementary  energy 
expression  and  the  unknown  constants  are  evaluated  by  minimising  the 


(A-9) 


At  =  an  arbitrary  small  change  associated  with  a  cariation  from  the 

co^Ltns'  specified  bo^ 

conditions  and  equilibrium  conditions  are  still  satisfied,  ^ 

V^*=  complementary  strain  energy  per  unit  volume,  and 

i  '  «“fi,r'  displacements 

The  expression  in  the  parentheses  is  defined  as  the  complementary  energy, 
tiiat  is, 


( A-  1C 
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complemeatary  energy.  The  stresses  obtained  will  approximately 
satisfy  the  compatibility  equations.  The  accuracy  of  the  approximation 
should  increase  with  the  addition  of  more  constants.  Since  the  compat¬ 
ibility  equations  are  not,  in  general,  exactly  satisfied,  the  accuracy  of 
any  displacements  subsequently  determined  from  such  a  solution  is  likely 
to  be  poor. 

The  Rayleigh- Ritz  method  is  applied  to  the  principle  of  minimum 
potential  energy  by  specifying  the  displacements  in  an  equation  form 
involving  undetermined  constants  such  that  the  boundary  conditions  given 
in  terms  of  displacements  are  satisfied.  Continuity  of  the  assumed  dis¬ 
placements  guarantees  that  the  compatibility  equations  will  be  satisfied. 

These  displacements  are  then  substituted  into  the  potential  energy 
expression  and  the  constants  determined  by  minimizing  the  potential 
energy.  The  displacements  obtained  by  substituting  the  values  of  these 
constants  into  the  assumed  functions  will  not,  in  general,  exactly  satisfy 
the  equilibrium  equations.  Even  though  this  method  can  provide  a  good 
approximation  to  the  displacements,  the  accuracy  of  the  associated 
stresses  is  likely  to  be  poor. 

In  either  case,  the  appropriate  energy  function  is  minimized  by 
equating  to  zero  the  derivative  of  the  energy  expression  with  respect  to 
each  unknown  constant,  that  is, 

=  0  ,  (A- 11) 

n 

where  E  is  the  energy  and  the  s  are  the  unknown  constants.  Such  a 
procedure  will  yield  n  equations  for  the  evaluation  of  the  c^*  s.  Normally, 
the  number  of  constants  is  increased  until  the  energy  approaches  a  sta¬ 
tionary  value.  When  this  point  is  reached,  the  addition  of  more  constants 
will  not  significantly  improve  the  approximation;  hence,  the  solution  is 
the  best  obtainable  with  the  chosen  approximations. 

It  should  be  noted  that  the  accuracy  of  the  solution  is  limited  by  the 
ability  of  the  assumed  stress  or  displacement  forms  to  approximate  the 
true  solution. 


Radial  deflection  of  inner  propellant  surface  for  112 propellant. 


Axial  displacement,  w,  in. 
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Fig.  4  Axial  deflec+lon  q]  end  (z*l)  of  propel lonf  groin  for  il2bw  propelbni 


Initial  distribution  of  this  report  was  made 
in  accordance  with  the  Joint  Army-  Navy- 
Air  Force  mailing  lists  for  Solid  Propellant 
and  Liquid  Propellant  technical  information 
plus  approved  supplements 


UNCLASSIFIED 


UNCLASSIFIED 


